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We study charged lepton flavour violation in low-scale seesaw models of minimal supergravity, which 
realize large neutrino Yukawa couplings thanks to approximate lepton-number symmetries. There 
are two dominant sources of lepton flavour violation in such models. The first source originates from 
the usual soft supersymmetry-breaking sector, whilst the second one is entirely supersymmetric 
>— ^ . and comes from the supersymmetric neutrino Yukawa sector. Within the framework of minimal 

supergravity, we consider both sources of lepton flavour violation, soft and supersymmetric, and 
calculate a number of possible lepton-flavour-violating transitions, such as the photonic decays of 
D , muons and taus, /i — ¥ &y, r — ► ey and r — >• fj,-y, their neutrinoless three-body decays, fi — > eee, 

pl«l 1 r —> eee, r — s- /i/i/i, r — > ee/i and r — ¥ e/i/i, and the coherent fx — ¥ e conversion in nuclei. 

After taking into account the exclusion bounds placed by present experiments of lepton flavour 
violation, we derive combined theoretical limits on the universal heavy Majorana mass scale mjv 
and the light-to-heavy neutrino mixings. Supersymmetric low-scale seesaw models offer distinct 
correlated predictions for lepton-flavour-violating signatures, which might be discovered in current 
and projected experiments, such as MEG, COMET/PRISM, Mu2e, super-BELLE and LHCb. 
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I. INTRODUCTION 



> 

On i 

Neutrino oscillation experiments [3-|3| have provided undisputed evidence of Lepton Flavour Violation (LFV) in 
the neutrino sector, pointing towards physics beyond the Standard Model (SM). Recent reactor neutrino oscillation 
experiments [|| have shown that the angle #13 of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix @ 
(N . is nonzero, thus hinting at a non-trivial neutrino-flavour structure and possibly at CP violation. Nevertheless, in spite 
of intense experimental searches 04H|, no evidence of LFV has been found yet in the charged lepton sector of the 
SM, implying conservation of the individual lepton numbers associated with the electron e, the muon fi and the tau r. 
All past and current experiments were only able to report upper limits on observables of Charged Lepton Flavour 
Violation (CLFV). 

Recently, the MEG collaboration [7| has announced an improved upper limit on the branching ratio of the CLFV 
decay /1 — > e^f, with B(n — > ej) < 2.4 x 10 -12 at the 90% confidence level (CL). As also shown in Table |H future 
experiments searching for the CLFV processes, /i — > ej, /i — > eee, coherent /1 — > e conversion in nuclei, r — > ej/n'-y, 
t — > 3 leptons and r — > lepton+light meson, are expected to reach branchin g-ra tio sensitivities to the level of 10 -13 [l2[ 
(If)" 14 Q), 10" 16 [3 (KT 17 Hg), 10" 17 llllii (KT 18 [SCLl), 10 -9 BlU, 10~ 10 and 10" 10 [3, respectively. 
The values in parentheses indicate the sensitivities that are expected to be achieved by the new generation CLFV 
experiments in the next decade. Most interestingly, the projected sensitivity for /i — > eee and \i — > e conversion in 
nuclei is expected to increase by five and six orders of magnitude, respectively. Given that CLFV is forbidden in the 
SM, its observation would constitute a clear signature for New Physics, thus rendering this field of investigations even 
more exciting. 

Although forbidden in the SM, CLFV is a generic feature for most of its extensions. One such well-motivated 
extension is the so-called Minimal Supersymmetric Standard Model (MSSM) [20[, where supersymmetry (SUSY) is 
softly broken at the 1-10 TeV scale for phenomenological reasons. The MSSM provides a quantum-mechanically stable 
solution to the gauge hierarchy problem and predicts rather accurate unification of the SM gauge couplings close to 
the Grand Unified Theory (GUT) scale. If R-parity is conserved, the lightest supersymmetric particle (LSP) is stable 
and, if neutral, such as the neutralino, it could play the role of the Dark Matter (DM) in the Universe. Finally, the 
MSSM typically predicts a SM-like Higgs boson lighter than 135 GeV, in agreement with the recent observations for 
a 125-GeV Higgs boson, made by the ATLAS [2l| and CMS [H| collaborations. 
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No. 


Observable 


Upper Limit 


Future Sensitivity 




1. 


B(p- 


■¥ erf) 


2.4 x 10" 12 [7] 


1 - 2 x 10 -13 [12], 10" 14 


[13] 


2. 


BQi- 


-> eee) 


10" 12 [8] 


lO" 16 [14], lO" 17 [13] 




3. 


pTi 




4.3 x 10" 12 [9], 


3 - 7 x 10" 17 [15, .16], 10 


- 18 [13, L7] 


4. 


dAu 




7 x 10" 13 [10] 


3 - 7 x 10" 17 [15, .16], 10 


" 1S [13.12] 


5. 


S(r - 


•>■ e 7 ) 


3.3 x 10" 8 [11] 


1 - 2 x 10" 9 [18, .19] 




6. 


B(r - 


-> |U7) 


4.4 x 10" 8 [11] 


2 x 10" 9 [18, .19] 




7. 


B{t - 


-> eee) 


2.7 x 10" 8 [11] 


2 x 10" 10 [18, .19] 




8. 


B(t - 


->• e/U/i) 


2.7 x 10" 8 [11] 


10" 10 [18] 




9. 


B{t - 


+ AWi) 


2.1 x 10 -8 [11] 


2 x 10 -10 [18,^9] 




10. 


B(r - 


> /zee) 


1.8 x 10" 8 [11] 


10" 10 [18] 





TABLE I: Current upper limits and future sensitivities of CLFV observables under study. 



In the MSSM with R-parity conservation, the lepton number is preserved and all left-handed light neutrinos v e ,n,T 
remain massless, exactly as in the SM. To account for the observed light neutrino masses and mixings, while maintain- 
ing R-parity, the field content of the MSSM needs to be extended. An interesting extension to the MSSM is provided 



by the so-called seesaw mechanism. There are three realizations of the seesaw mechanism: the seesaw type I [23[ , the 
seesaw type II [24[ and the seesaw type III [25[ . The three scenarios differ by the nature of their seesaw messengers 
that are needed to explain the small neutrino masses. In this study, we will adopt a low-scale variant of the seesaw 
type-I realization, whose messengers are three singlet heavy neutrinos Ni t 2,z- 

In the usual seesaw type-I mechanism, the heavy singlet neutrinos must assume masses of order ~ 10 12 ~ 14 GeV, for 
electroweak-scale Dirac neutrino masses, in order to account for the observed light neutrino mass spectrum. The mixing 
between light and heavy neutrinos is of the order £„jv ~ \J m v /mN ~ 10~ 12 , for light-neutrino masses m„ ~ 10 _1 eV. 
As a consequence, heavy neutrinos decouple from low-energy processes of CLFV in the SM with right-handed neutrinos, 
giving rise to extremely suppressed and unobservable rates. In the MSSM with right-handed neutrinos, however, the 
singlet heavy neutrinos do not fully decouple. They impact the low-energy sector, through renormalization group (RG) 
effects that induce sizeable LFV in the slepton sector for a SUSY mass scale Msusy ~ 1 — 10 TeV. 

A potentially interesting alternative to the ordinary seesaw mechanism may arise from the presence of approximate 
lepton-number symmetries [26l - |29j in the theory. The smallness of the light neutrino masses is a consequence of these 
approximate leptonic symmetries which are radiatively stable [29l . |30j , whilst the heavy neutrino mass scale could be 
as low as 100 GeV. Unlike in the usual seesaw scenario, the light-to-heavy neutrino mixings £„jy are not correlated 
to the light neutrino masses m v . Instead, £„jv are free parameters, constrained by experimental limits on deviations 
of the W - and Z-boson couplings to leptons with respect to their SM values [3ll l32j j . Approximate lepton-number 
symmetries do not restrict the size of LFV, and so potentially large phenomena of CLFV may be predicted. This 
feature is quite generic both in the SM [33| and in the MSSM [lj,[3^] augmented with low-scale right-handed neutrinos. 
It is this new source for LFV in the MSSM that we wish to study in detail here, in addition to the one resulting from 
the frequently considered soft SUSY breaking sector [36| - |39| . 

In this article, we denote for brevity the SM and the MSSM extended by low-scale right-handed neutrinos and 
approximate lepton-number symmetries by ^rSM and j'rMSSM, respectively. Our study in this paper focuses on 
the z/rMSSM with constrained boundary conditions at the gauge-coupling unification scale, within the framework 
of minimal supergravity (mSUGRA). However, the results presented here are applicable to more general soft SUSY- 
breaking scenarios. 

The z/rMSSM has some interesting features with respect to the MSSM. In particular, the heavy singlet sneutrinos 
may emerge as new viable candidates of Cold Dark Matter (CDM) [40J. In addition, the mechanism of low-scale 
resonant leptogenesis [30, |4l| could provide a possible explanation for the observed Baryon Asymmetry in the Uni- 
verse (BAU), as the parameter space for successful electroweak baryogenesis gets squeezed by the current LHC 
data |42| . 

Given the multitude of quantum states mediating LFV in the j^rMSSM, the predicted values for observables of 
CLFV in this model turn out to be generically larger than the corresponding ones in the i<rSM, except possibly for 
B(l — > I'j) [UHII, where I, V — e, /Lt, r. The origin of suppression for the latter branching ratios may partially 
be attributed to the SUSY no-go theorem due to Ferrara and Remiddi (43|, which states that the magnetic dipolc 
moment operator necessarily violates SUSY and it has therefore to vanish in the supersymmetric limit of the theory. 

The goal of this paper is to improve upon earlier calculations of CLFV in supersymmetric low-scale seesaw models 
which are either approximate, or incomplete. Specifically, in [34| a low-scale aligned neutrino and sneutrino spectrum 



3 



was assumed to the leading order in the heavy neutrino mass, whilst in [44|, |45[ only a subset of loop effects mediated 
by the photon 7, the Z boson, the neutral Higgs scalars and box graphs was considered. Instead, we present here 
detailed analytical expressions for CLFV observables, induced by 7- and Z-boson formfactors, as well as the complete 
set of box contributions involving heavy neutrinos, sleptons, sneutrinos, charginos, neutralinos and charged Higgs 
bosons. We also derive the complete set of kinematic formfactors that contributes to the three-body LFV decays of 
the muon and tau, such as /1 — > eee, r — » eee, r — > r — > ee/i and r — > e\x\x. In our numerical analysis, we 

consider benchmark scenarios of the i/rMSSM, which are in agreement with the existing LHC data from the recent 
discovery of a SM-like Higgs boson (21,|22| and the non-observation of squarks and gluinos [46j . 

The paper is organized as follows. In Section |H] we describe the relevant leptonic sector of the ^_rMSSM and 
introduce two baseline scenarios based on approximate lepton- number symmetries. Section [IIII contains details of the 
calculation of the LFV amplitudes and branching ratios of CLFV decays of the tau and muon, where the complete 
set of chiral structures of the amplitudes contributing to the CLFV processes / — > Z'7, I — > l'l\li and /i — > e conversion 
in nuclei are derived. Section [IVI presents numerical estimates for the aforementioned processes of CLFV, within the 
two baseline mSUGRA scenarios introduced in Section [TTJ Moreover, we discuss correlated predictions of the CLFV 
observables with B(l — > i'7) and other kinematic parameters. Section IVl summarizes the results of our analysis and 
presents our conclusions. All technical details have been relegated to the appendices. In Appendix [5] we give all 
relevant interaction vertices in the i^jMSSM. Appendix [B] defines the one-loop functions that appear in our analytic 
calculations. In terms of these loop functions, Appendix[C]describes the analytic results for all the one- loop formfactors 
that occur in the CLFV transition amplitudes. 



II. LOW-SCALE SEESAW MODELS AND LEPTON FLAVOUR VIOLATION 



In this section we describe the leptonic sector of the MSSM extended by low-scale right-handed neutrinos, which 
we call the z/rMSSM for brevity. In addition, we introduce the neutrino Yukawa structure of two baseline scenarios 
based on approximate lepton-number symmetries and universal Majorana masses at the GUT scale. We will use these 
two scenarios to present generic predictions of CLFV within the framework of mSUGRA. 

The leptonic superpotential part of the i^rMSSM reads: 

Wiepton = E c h e H d L + N c h u LH u + ^N c m M N c , (2.1) 

where H. u ,di L, E and N c denote the two Higgs-doublet superfields, the three left- and right-handed charged- lepton 
superfields and the three right-handed neutrino superfields, respectively. The Yukawa couplings h CyL , and the Majorana 
mass parameters tocim form 3x3 complex matrices. Here, the Majorana mass matrix raj/ is taken to be SO (3)- 
symmetric at the GUT scale, i.e. mj^ = 1 3 . 

In this study, we consider low-scale seesaw models p6l - l29| . where the smallness of the light neutrino masses is 
protected by natural, quantum-mechanically stable cancellations due to the presence of approximate leptonic symme- 
tries [29], whilst the Majorana mass scale tun can be as low as 100 GeV. In these models, the neutrino induced LFV 
transitions from a charged lepton I — fi, t to another charged lepton V ^ I are functions of the ratios [H, H3, H§[ 

2 3 
= ^\{hlK) in = y2B VNi B lNi (2.2) 

Zm N i=1 

and are not constrained by the usual seesaw factor m v /mN : where v u /^/2 = (H u ) is the vacuum expectation value 
(VEV) of the Higgs doublet H u , with tan/3 = (H u ) j (Hj) . Moreover, the mixing matrix Bi^ i that occurs in the 
interaction of the bosons with the charged leptons / = e, [i, r and the three heavy neutrinos Vi ; 2,3 is defined in 
Appendix |A"1 Note that the LFV parameters fli'i do not directly depend on the RG evolution of the soft SUSY-breaking 
parameters, except through the VEV v u at the minimum of the Higgs potential. 

In the weak basis {{v e .n.TL) C , ^1.2,3^}, the neutrino mass matrix in the ^rMSSM takes on the standard seesaw 
type-I form [23] : 

M„ = f ° T m f ) , (2.3) 

where = sin/? g^hl and are the Dirac- and Majorana-neutrino mass matrices, respectively. In this 

paper, we consider two baseline scenarios of neutrino Yukawa couplings. The first one realizes a U(l) leptonic 
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symmetry [30| and is given by 



to 





ae~~ be~ 





ce^ 



(2.4) 



In the second scenario, the structure of the neutrino Yukawa matrix h„ is motivated by the discrete symmetry group 
A4 and has the following form (4t|: 




(2.5) 



In (|2.4[) and (|2.5[) . the Yukawa parameters a, 6 and c are assumed to be real. If the above leptonic (discrete or 
continuous) symmetries are not broken, the light neutrinos are massless. The neutrino masses are obtained by adding 
small terms that break the symmetry of the Yukawa matrix. It is essential to remark here that the predictions for 
CLFV are independent of the flavour structure of these small symmetry-breaking terms in our low-scale seesaw models 
under study. These terms are needed to fit the low-energy neutrino data. For this reason, we do not discuss here 
particular symmetry breaking patterns of the above two baseline Yukawa scenarios given in (|2.4I) and (12.51) . 

The second source of LFV in the models under consideration originates from scalar- neutrino (sneutrino) interactions, 
namely the supersymmetric partners of the left-handed and right-handed neutrinos. Specifically, the sneutrino mass 
Lagrangian in flavour and mass bases is given by 



( * \ 



cF = 



/~t ~G\ 



~T ~CT 



)M% 



(2.6) 



where M? is a 12 x 12 Hermitian mass matrix in the flavour basis and M| is the corresponding diagonal mass matrix 
in the mass basis. More explicitly, in the flavour basis {v e ,n,TL 
M§ may be cast into the form: 



rL'^i2 3fl}' the sneutrino mass matrix 



M? 



/ Hi 
Nt 


V Mt 

where the block entries are the following 3x3 matrices 



N 

H 2 

M* 






M T 



M \ 


N* 
H 2 J 



(2.7) 



Hi = rri£ + mDm' D + -M|cos2^ 



H 2 
M 
N 



,t 



m s + m' D m D + m M m M 
m £) (A, y - /icot/3) 
m D m M . 



(2.8) 



Here m|, m? and A„ are 3x3 soft SUSY-breaking matrices associated with the left-handed slepton doublets, the 
right-handed sneutrinos and their trilinear couplings, respectively. In the supersymmetric limit, all the soft SUSY- 
breaking matrices are equal to zero, tan j3 = 1 and /i = 0. As a consequence, the sneutrino mass matrix M? can be 
expressed in terms of the neutrino mass matrix M„ in (|2.3I) as follows: 



Mi 



SUSY 



M„Mt 6x6 
6x6 MtM„ 



(2.9) 



leading to the expected equality between neutrino and sneutrino mixings. Sneutrino LFV mixings do depend on the 
RG evolution of the i/rMSSM parameters, but unlike the LFV mixings induced by soft SUSY-breaking terms, the 
sneutrino LFV mixings do not vanish at the GUT scale. 
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FIG. 1: Feynman graphs contributing to / — > I'j and Z I /'(£—» ZZ') amplitudes 



The sneutrino LFV mixings are obtained as combinations of unitary matrices that diagonalize the sneutrino, slepton 
and chargino mass matrices. It is interesting to notice that in the diagonalization of the sneutrino mass matrix M? 
in (12. 7[) . the sneutrino fields £ e ,/i,TL) £p 2 3R and their complex conjugates £* Mr £, 3i? are treated independently. 
As a result, the expressions for £ e ,^,TL and z/p 2 3fl , in terms of the real-valued mass eigenstates -/V lj2 127 arc n °t 
manifestly complex conjugates to v* L and v^zr: thus leading to a two-fold interpretation of the flavour basis 
fields, 

v* = &)* = wfiiv^ , 

= Wf +6 ^ , (2.10) 

where Pi 2,3 = v e ,\x,TL an d £4,5.6 = ^p2 3i?> w ith % = 1, 2, . . . , 6 and A = 1, 2, ... , 12. For this reason, in Appendix |A1 
we include all equivalent forms that the Lagrangians, such as and Cjjjqz, can be written down. 

One technical comment is now in order. Unlike [13, 111], the diagonalization of the 12 x 12 sneutrino mass ma- 
trix M| and the resulting interaction vertices are evaluated numerically, without approximations. To perform the 
diagonalization of M? numerically, we use the method developed in (Hoj for the neutrino mass matrix. This method 
becomes very efficient, if the diagonal submatrices have eigenvalues larger than the eigenvalues of the other subma- 
trices. Therefore, we assume that the heavy neutrino mass scale tun is of the order of, or larger than the scale of the 
other mass parameters in the ^MSSM. 

Finally, a third source of LFV in the i/rMSSM results from soft SUSY-breaking LFV terms. These LFV terms are 
induced by RG running and are absent at the GUT scale in mSUGRA. Their size strongly depends on the interval of 
the RG evolution from the GUT scale to the universal heavy neutrino mass scale to at. 

All the three different mechanisms of LFV, mediated by heavy neutrinos, heavy sneutrinos and soft SUSY-breaking 
terms, depend explicitly on the neutrino Yukawa matrix h„ and vanish in the limit h, y — > 0. 
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III. CHARGED LEPTON FLAVOUR VIOLATION 

In this section, we present key details of our calculations for a number of CLFV observables in the ^_rMSSM. In 
detail, Section fill Al gives analytical results for the amplitudes of the CLFV decays: I — > I'j and Z — > ll' c , and their 
branching ratios. Correspondingly, Sections MI Bl and MI CI give analytical expressions for the neutrinoless three-body 
decays I — > l'lil§ pertinent to muon and tau decays, and for coherent /i — > e conversion in nuclei. All these analytical 
results are expressed in terms of one-loop functions and composite formfactors that are defined in Appendices [Bl 
andO 

A. The Decays I l'f and Z ll'° 

At the one-loop level, the effective couplings jl'l and Zl'l are generated by the Feynman graphs shown in Fig. [1] 
The general form of the transition amplitudes associated with these effective couplings is given by 



(F^)m (q 2 ^ - iq^)P L + (F?) vl (q% - iq»)P h 



+ {G^)vii°^q U PL + {GK)i,iio^q v PR I, (3.1) 



q-Zl'l _ 9W <XW J, 

M 8n cos 9 W 



{F%) vlllx P L + (Fi) inltl P R I, (3.2) 



where Pl(r) = 5 [1 — (+)7s]> a w = fl^/G^Tr), e is the electromagnetic coupling constant, Mw = gwy/v 2 + v d/^ 1S tne 
MK-boson mass, 6 W is the weak mixing angle and q = pi> — pi is the photon momentum. The form factors (F^)i'i, 

(F^)i'i (G^)i'i, (G~)i>i, (Fz)i>i and (Fz)i'i receive contributions from heavy neutrinos Ni^s, heavy sneutrinos Ni >2 ,3 
and RG induced soft SUSY-breaking terms. The analytical expressions for these three individual contributions are 
given in Appendix O Note that according to our normalization, the composite formfactors (G^)vi and (G^)in have 
dimensions of mass, whilst all other formfactors are dimensionless. 

It is essential to remark here that the transition amplitudes (|3. II) and (|3.2|) are also constituent parts of the leptonic 
amplitudes / — > l'lil§ and semileptonic amplitudes I — > l'qiq~2, which will be discussed in more detail in Sections MI Bl 
and MI Cl To calculate the CLFV decay I — > I'j, we only need to consider the dipole moment operators associated 
with the formfactors (G^)in and (G^)in in (|3.1I) . Taking this last fact into account, the branching ratios for I —> l >r y 
and Z -> ll' c + l c V are given by 

B ^^ = #]i^(l^'| 9 + l( G >'| a )' (3 - 3) 
B{Z^l'l c + l' c l) = -^^-(\ { F L z ) vl \ 2 + \{Fi) vl \ 2 ) . (3.4) 

Observe that the above expressions are valid to leading order in external charged lepton masses and external momenta, 
which constitutes an excellent approximation for our purposes. Thus, in (I3.4[) we have assumed that the Z-boson 
mass Mz is much smaller than the SUSY and heavy neutrino mass scales, Msusy and tojv, and we have kept the 
leading term in an expansion of small momenta and masses for the external particles. In the decoupling regime of 
all soft SUSY-breaking and charged Higgs-boson masses, the low-energy sector of the z^rMSSM becomes the VrSM. 
In this z/rSM limit of the theory, the analytical expressions for B(l — > Z'7) and B(Z — > l'l c + l' c l) take on the forms 
presented in [5l| and [47|, respectively. 



B. Three-Body Leptonic Decays I — >• l'l±l 



We now study the three-body CLFV decays / — > l'lil§ > where I can be the muon or tau, and I', l\, I2 denote other 
charged leptons to which I is allowed to decay kinematically. 

The transition amplitude for I — > I'lxl^ receives contributions from 7- and Z-boson-mediated graphs shown in Fig.[T] 
and from box graphs displayed in Fig. [21 The amplitudes for these three contributions are 
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FIG. 2: Feynman graphs contributing to the box / — > Z'ZiZ 2 amplitudes. 
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where g l L = — 1/2 + and = are Z-boson-lepton couplings and s w = sin^. The composite box formfactors 
Bf A Y , A = V, S, T, X, Y = L, R are given in Appendix [Cl The labels V, S and T denote the formfactors of the vector, 
scalar and tensor combinations of the currents, while L and R distinguish between left and right chiralities of these 
currents. We note that the box formfactors contain both direct and Fierz-transformed contributions. Equation (|3.8[) 

represents a shorthand expression that takes account of all individual contributions to the amplitude T^ox 1 2 induced 



by box graphs. Explicitly, the matrices appearing m 



read: 



t^L t^R t^L t^R j^L j^H 
1 V: 1 V: 1 S> 1 St 1 T> 1 T 



(lfj.PL, JuPr, Pl, Pr, v^vPl, o~hvPr 



(3.9) 



Notice that the tensor formfactors -Bfc? and -BSf vanish in the sum 



i.e. -B^? = Bp^ = 0, as a consequence of 
pT . A very similar chiral structure holds also true for the semileptonic box amplitudes 



B RL 

IT 



do not include in their calculations the chiral 



the identity 0^75 

defined in the next section. We should remark here that [37jj and 
structures Pl x Pr and Pr x Pl and their corresponding formfactors B\g and Bfg. 

In a three-generation model, the transition amplitude for the decays I — > I'lxl^ may fall in one of the following three 
classes or categories [Hj]: (i) V ^ l\ = I2, (ii) I' = h = h, and (iii) V = h 7^ li. In the first two classes (i) and (ii), 
the total lepton number is conserved, whilst in the third class (iii) the total lepton number is violated by two units. 
Here, we ignore this lepton-number violating class (iii) of charged lepton decays, since the predictions turn out to be 
unobservably small in the ^rMSSM. Moreover, we suppress the universal indices VI that appear in the photon and 
Z-boson formfactors, i.e. F^, F R , F^ and F%. Given the above simplification and the notation of the box formfactors 
in (|3.8I) . the branching ratios for the class (i) and (ii) of CLFV three-body decays are given by 



B(l -> I'hlf) = 



"l "to 
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(3.10) 
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(I' 




1 R-R-R12 



(3.11) 



where m; and me, m; 2 are the masses of the initial- and final-state charged leptons and T/ is the decay width of 
the charged lepton I. Here we should emphasize that the transition amplitudes (|3.5p . (|3.6p . p.7| and the branching 
ratios (|3.10p and (|3.1ip have the most general chiral and formfactor structure to leading order in external masses and 
momenta and so they are applicable to most models of New Physics with CLFV. Finally, we have checked that the 
branching ratios (|3.10p and (|3.11|> go over to the results of [33], in the i/rSM limit of the theory. 
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C. Coherent /x — 5- e Conversion in a Nucleus 



The coherent /i — > e conversion in a nucleus corresponds to the process J M — > e~J + , where J^ t is an atom of 

nucleus J with one orbital electron replaced by a muon and J + is the corresponding ion without the muon. The 
transition amplitude for such a CLFV process, 

T" eJ = (J + e-\T d ^ de \J„) + (J + e-\T u ^ ue \J„} , (3.12) 

depends on two effective box operators, 

r b d r de = -Jr E E ^/er^Jrid 

W X,y=L,fl A=V,S,T 

_ 2 

(<^d) e (V d R P R + V d L P L )n, (3.13) 



2M^ 



'box - TVp- 2^ 2^ ^ eL A nul A u 



1 II* 



2M^ 



(u+u) e (K* P fl + K L Pl) ■ (3-14) 



Here /j, and e are the muon and electron wave functions and d and u are field operators acting on the and J + 
states. The formfactors B dA Y and B* A are given in Appendix ICl The composite formfactors V^ L , V^, V R , V R may 
conveniently be expressed as 



VJ- = -l^-^) + (\-\^)Fk + \(B^ + B^ + B^ + B^ 



V U L = \sl (l% - ±G R ) + ( - \ + H s 2 ) F| + i + itf * + B*? + B«# 



Kf = §4 - — G^j + ( - - + |4 J ^ + j ( + + B«s + S«s ) , (3-15) 



where , F R , F%, F§ is the shorthand notation for (F^) ell , (F R ) ef ^ (F|) eM! (F R ) efl . 

Our next step is to determine the nucleon matrix elements of the operators vJu and cftd. These are given by 

(J+e-|utu|J„) = (2Z + N)F(-ml) , 

(J + e-\dU\J„) = (Z + 2N)F{-ml) , (3.16) 

where the formfactor F(q 2 ) incorporates the recoil of the J + ion (53|, and the factors 2Z + N and Z + 2N count the 
number of u and d quarks in the nucleus J, respectively. Hence, the matrix element for — > J + /i~ can be written 
down as 



with 



T J ^ J e = -^TF(-ml)e(QwPR + QwPL)v, (3.17) 



Qw = (2Z + N)Vf + {Z + 2N)V d L 



)« - (2Z + N)V* + (Z + 2N)V d R . (3.18) 



Given the transition amplitude (|3.17l) . the decay rate J M — » J + e is found to be 



" e 167r 2 M^r capture Z 



\n-^l)\ 2 {\Qw\ 2 + \QU) . ( 3 - 19 ) 



where r captU rc is the capture rate of the muon by the nucleus, and Z e g is the effective charge which takes into account 
coherent effects that can occur in the nucleus J due to its finite size. In our analysis, we use the values of Z e s quoted 
in 54| . We reiterate that the branching ratio (13. 19)) possesses the most general formfactor structure to leading order 



in external masses and momenta and is relevant to most models of New Physics with CLFV. Finally, we have verified 
that our analytical results are consistent with [34| [5^, [56[ in the iVrSM limit of the theory. 
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IV. NUMERICAL RESULTS 

In this section, we present a numerical analysis of CLFV observables in the i/rMSSM. In order to reduce the number 
of independent parameters, we adopt the constrained framework of mSUGRA. In detail, our model parameters are: 
(i) the usual SM parameters, such as gauge coupling constants, the quark and charged-lepton Yukawa matrices 
inputted at the scale Mz, (h) the heavy neutrino mass to at and the neutrino Yukawa matrix h„ evaluated at m^r, 
(iii) the universal mSUGRA parameters m , M 1 / 2 and A inputted at the GUT scale, and (iv) the ratio tan/3 of the 
Higgs VEVs and the sign of the superpotential Higgs-mixing parameter /i. 

The allowed ranges of the soft SUSY- breaking parameters mo, M1/2, Aq and tan/3 are strongly constrained by a 
number of accelerator and cosmological data [U [52, 5(| [57} • For definiteness, we consider the following set of input 
parameters: 

tan /3= 10, mo = 1000 GeV , A = -3000 GeV , M 1/2 = 1000 GeV . (4.1) 

Here we take the /x parameter to be positive, whilst its absolute value \fi\ is derived form the minimization of the 
Higgs potential at the scale Mz- With aid of [58l - l60l |. we verify that the parameter set ( 14.11) predicts a SM-like Higgs 
boson with mu ~ 125 GeV, in agreement with the recent discovery at the LHC [U |22|. and is compatible with the 
current lower limits on gluino and squark masses (46[. The set (|4.1I) is also in agreement with having the lightest 
neutralino as the Dark Matter in the Universe [57j . 

We employ the one-loop RG equations of [HI, [62| to evolve the the gauge coupling constants and the quark and 
charged lepton Yukawa matrices from Mz to the GUT scale, while the heavy neutrino mass matrix vclm and the 
neutrino Yukawa matrix h. v are evolved from the heavy neutrino mass threshold to at to the GUT scale. Furthermore, 
we assume that the heavy neutrino-sneutrino sector is supersymmetric above to at. For purposes of RG evolution, 
this is a good approximation for to at larger than the typical soft SUSY-breaking scale [34|. At the GUT scale, the 
mSUGRA universality conditions are used to express the soft SUSY-breaking masses, in terms of too, Mi/ 2 and A$. 
Hence, all scalar masses receive a soft SUSY-breaking mass mo, all gaugino are mass-degenerate to Mi/ 2 , and all scalar 
trilinear couplings are of the form \i x Aq, with x — u,d, I, v, where are the Yukawa matrices at the GUT scale. The 
sneutrino mass matrix acquires additional contributions from the heavy neutrino mass matrix. The sparticle mass 
matrices and trilinear couplings are evolved from the GUT scale to Mz, except for the sneutrino masses which are 
evolved to the heavy neutrino threshold to at. Having thus obtained all sparticle and sneutrino mass matrices, we can 
numerically evaluate all particle masses and interaction vertices in the i^rMSSM, without approximations. 

To simplify our numerical analysis, we consider two representative scenarios of Yukawa textures discussed in Sec- 
tion [H] Specifically, the first scenario realizes the U(l)-symmetric Yukawa texture in (|2.4p . for which we take either 
a = b and c = 0, or a — c and b — 0, or b — c and a — 0, thus giving rise to CLFV processes /1 — > eX, t — > eX 
and r — > /j,X, respectively. Here X represents the lepton-flavour conserving state(s), e.g. X = 7, e + e~, /x + fj," , qq. 
The second scenario is motivated by the A4 group and uses the Yukawa texture (|2.5p . where the parameters a, b and c 
are taken to be all equal, i.e. a = b = c. 

The heavy neutrino mass scale toat strongly depends on the size of the symmetry-breaking terms in the Yukawa 
matrix h„. For instance, for the model (|2.4p . the typical values of the [/(l)-lepton-symmetry-breaking parameters 
ei = t e ,ii,T consistent with low-scale resonant leptogenesis is e < 10 -5 [30| . leading to light-neutrino masses 

m „ „ if. „ lO-eV^)' (iM) . ,4,, 
tojv MO b / V rriN / 

Taking into account the constraint m v 10 _1 eV generically derived from neutrino oscillation data, we may estimate 
that the heavy neutrino mass scale to^t is typically restricted to be less than 10 TeV, for e; = 10~ 5 . If the assumption 
of successful low-scale leptogenesis is relaxed, the symmetry-breaking parameters e; has only to be couple of orders 
in magnitude smaller than the Yukawa parameters a, b and c, with a, b, c ^ 10. Thus, for e; < 10~ 3 — 10~ 2 , the 
heavy neutrino mass scale ttin may be as large as 10 7 — 10 9 TeV, leading to the decoupling of heavy neutrinos from 
low-energy observables. As our interest is in the interplay between heavy neutrino, sneutrino and soft SUSY-breaking 
contributions to CLFV observables, we will only study here the parameter space in which to^ < 10 TeV. 

In the present analysis, we consider that the symmetry-preserving Yukawa parameters a, b and c are limited through 
the perturbativity condition: Tr hj/ h. v < An, which we require to hold true for the entire interval of the RG evolution: 
\n(Mz/TcV) < t < ln(MGUT/TeV) . For the model in (|2.4[) , this condition translates into the constraint: a < 0.34, 
and for the model in (|2.5p . to a < 0.23. As a consequence, we do not display in plots numerical values for points in 
parameter space, for which the aforementioned perturbativity condition gets violated. 

In Fig. [21 we display numerical predictions for the fi-LFV observables B(/i — > eX): B(/i —> ej) [blue (solid) line], 
B([i —> eee) [red (dashed) line], R^ e [violet (dotted) line] and [green (dash-dotted) line], as functions of B(pL — > ej) 
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B(yU->ey) a=b,c=0 




B(yU->ey) a=b=c 



FIG. 3: Numerical estimates of B(p —¥ ey) [blue (solid)], B(ft eee) [red (dashed)], R^l [violet (dotted)] and R^e [green 
(dash-dotted)], as functions of B(/i — > ej) (left pannels) and the Yukawa parameter a (right pannels), for mjv = 400 GeV and 
tan/3 = 10. The upper two pannels correspond to the Yukawa texture (|2.4[) . with a — b and c = 0, and the lower two pannels 
to the Yukawa texture ()2.5[) . with a — b — c. 



(left pannels) and the Yukawa parameter a (right pannels), for m^r = 400 GeV and tan j3 = 10. The upper two pannels 
assume the Yukawa texture in (|2.4[) . with a = b and c = 0, whilst the lower two pannels correspond to the Yukawa 
texture in (|2.5I) . with a = b = c. In Fig. [U we give numerical estimates for the same set of /.t-LFV observables, but 
for tun = 1 TcV. In Figs. [3] and |H the Yukawa parameter a has been chosen, such that 10 -20 < B(fi — > ey) < 10~ 10 . 
Such a range of values includes both the present 0413, IH, H3] and future [H, [13, EH Ir35rl68| experimental limits. 
As we see from Figs. [3] and [U the CLFV observables under study depend quadratically on the Yukawa parameter a, 
namely they are proportional to a 2 . Instead, the quartic Yukawa terms proportional to a 4 [34j remain always small, 
which is a consequence of the imposed perturbativity constraint: Tr(hJ,h„) < An, up to the GUT scale. 

By analogy, Figs. [5] and present numerical estimates of the t-LFV observables B(t — > eX): B(r — > erf) [blue 
(solid) lines], B(t — > eee) [red (dashed) lines] and B(t — > e\i[i) [violet (dotted) lines], as functions of B(t — > erf) 
(left pannels) and the Yukawa parameter a (right pannels), for tojv = 400 GeV and mjv = 1 TeV, respectively. Note 
that we do not show predictions for the fully complementary observables B(t — > /J.X): B(t — > fJ/y), B[t — > fififi) and 
B(t — > fiee). The upper pannels give our predictions for the Yukawa texture (|2.4[) . with a = c and 6 = 0, and the 
lower pannels for the Yukawa texture (|2.5p . with a = b = c. In both Figs. and the Yukawa parameter a has been 
chosen, such that 10~ 16 < B(r — > ey) < 10 -7 . As can be seen from Figs. [5] and O all observables B(t — > eX) of 
t-LFV (with X = 7, ee, /x/i) exhibit similar quadratic dependence on the small Yukawa parameter a. However, close 
to the largest perturbatively allowed values of a, i.e. a 5 0.34 for the model (|2.4[) and a < 0.23 for the model (|2.5[) . 
some of the observables of t-LFV exhibit either numerical instability, or the existence of a cancellation region in 
parameter space, as we will see below. 

Figure [7] presents numerical estimates of B(fi — > ey) [blue (solid) line], B(fi —> eee) [red (dashed) line], [violet 
(dotted) line] and [green (dash-dotted) line], as functions of B([i — > ey) (left pannels) and the heavy neutrino mass 
scale rriN (right pannels). In all pannels, we keep the Yukawa parameter a fixed by the condition B{p — > ey) = 10 -12 
for tow = 400 GeV, using the benchmark value: tan/? = 10. The upper pannels display numerical values for the 
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B(yU->ey) a=b,c=0 




B(yU->ey) a=b=c 



FIG. 4: The same as in Fig.O but for m N = 1 TeV. 

Yukawa texture (|2.4[) . with a = 6 and c = 0, and the lower pannels for the Yukawa texture (|2.5[) , with a = = c. 
The heavy neutrino mass is varied within the LHC explorable range: 400 GeV < win < 10 TeV. All observables 
.B(/z — ► eX) of /i-LFV (with X = 7, ee, Ti, Au) exhibit a non-trivial dependence on mjv- The branching ratio 
B(n — > ej) shows a dip at rnjv ~ 800 GeV in both models (|2.4I) and (|2.5I) . signifying the existence of a cancellation 
region in parameter space, due to loops involving heavy neutrino, sneutrino and soft SUSY-breaking terms. For 
win ~ 3 TeV, all observables tend to a constant value, as a result of the dominance of the soft SUSY-breaking 
contributions. 

In Fig. [3] we show contours of the Yukawa parameters (a, 6, c) versus the heavy neutrino mass scale rriN, for 
B([i — > ej) [blue (solid) line], B([i—> eee) [red (dashed) line], iff* [violet (dotted) line] and iZ„" [green (dash-dotted) 
line]. The Yukawa parameter a and mjv are determined by the condition B([i — > ej) = 10~ 12 . The labels in the 
vertical axes indicate the two Yukawa textures in (12.41) and (|2.5I) . which we have adopted in our analysis. The contours 
for B(fi — > ey) display a maximum for ijin ~ 800 GeV, as a consequence of cancellations between heavy neutrino, 
sneutrino and soft SUSY-breaking contributions (cf. Fig. [7]). 

Figure [3] shows contours of the Yukawa parameters (a,6,c), as functions of mjv, for B{t — > ej) [blue (solid) line], 
where the parameters a and m^r are determined by the condition B(t — > cy) = 10 -9 . We do not give numerical 
results for B(t — > ^7), as these are fully complementary to the ones given for B(t — > &y). Notice that given the above 
condition on B(t —¥ ej), no solution exists for B(t — > eee) and B(t — > efifi). 

In our numerical analysis so far, we have kept the value of tan/3 fixed to its benchmark value given in (|4.1I) : 
tan /3 = 10. In Fig. 1101 we relax this assumption, varying tan/3 in the interval 5 ^ tan/3 < 20, while maintaining 
agreement with a SM-like Higgs boson mass Mh ~ 125 GeV and taking into account that the combined experimental 
and theoretical errors are currently of the order of 5-6 GeV. Specifically, in Fig. [TU] we display the dependence of 
B([i — > ej) [blue (solid) line], B(fi — > eee) [red (dashed) line], R^l [violet (dotted) line] and [green (dash-dotted) 
line] on tan/3. In all pannels, the Yukawa parameter a is determined by the condition B{ji — > ej) = 10~ 12 . The upper 
pannels in Fig. [TOl show numerical results for m^y = 400 GeV, while the lower pannels for = 1 TeV. The left 
pannels give our predictions for the Yukawa texture (I2.4p . with a — b and c = 0, and the right pannels for the Yukawa 
texture (|2.5p . with a = b = c. In the lower pannels, we observe a suppression of B(/i —¥ ej), for values tan/3 « 7, due 
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B(r->ey) a=b=c 



FIG. 5: Numerical estimates of B(t — > ej) [blue (solid)], B(t — > eee) [red (dashed)] and B(t — > efifi) [violet (dotted)], as 
functions of B(r — > ey) (left pannels) and the Yukawa parameter a (right pannels), for rrejv = 400 GeV and tan/3 = 10. The 
upper pannels present predictions for the Yukawa texture (|2.4[) . with a = c and 6 = 0, and the lower pannels for the Yukawa 
texture (12. 5|) . with a = b — c. 



to cancellation between heavy neutrino, sneutrino and soft SUSY-breaking effects. 

It is interesting to compare the contributions of the magnetic dipole formfactors to the CLFV observables, with those 
originating from the remaining formfactors. Specifically, if one assumes that only the magnetic dipole formfactors G^' R 
contribute in (|3.10p . p. lip and p,19j) . then the following analytical results are obtained for the ratios: 

Bjl : -> VWlt) = ^4^) (4 .3) 
B(l -> 3tt V mf, J ' 

R 2 = "> m,C ) = -fifing -^ (44) 



B(Z Z'7) 



i? 3 = = 16a 4 -^ZZ e 4 ff |F(- M 2 )| 2 . (4.5) 

5(/i -> e7j r cap turo 

According to the formulae (|4.3p - (|4.5p . the predicted R\ values for t — > e/x/j, and r — > fj,ee are 1/90 and 1/419 
respectively, the predicted i?2 values for /i — > eee, r — > eee and the r — > are 1/159, 1/91 and 1/460 respectively, 
and the predicted i?3 values for Ti and Au are 1/198 and 1/188 respectively. 

In Fig. [TTJ we give numerical estimates of the ratios Ri{\i — > eee), Rj 1 and R% u , as functions of mjy. The Yukawa 
parameter a is fixed by the condition B(fi —> ey) — 10~ 12 , for tun — 400 GeV and tan/3 = 10. In the upper pannel, 
thick lines show the predicted values obtained by a complete evaluation of R^ijJL — > eee) [blue (solid) line], Rj l [red 
(dashed) line] and R^ u [violet (dotted) line] , while the respective thin lines are obtained by keeping only the magnetic 
dipole formfactors and G R . Hence, we see that going beyond the magnetic dipole moment approximation may 
enhance the ratios i?2,3 by more than two orders of magnitude. 

The two middle pannels of Fig. [TT] provide a formfactor analysis of i?2(^ eee) and R$ u , by considering separately 
the contributions due to G 1 and F 7 [blue (solid) line] , Fz [red (dashed) line] and box formfactors [violet (dotted) line] . 
In particular, we observe that heavy neutrino contributions to the box formfactors become comparable to and even 
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B(r->ey) a=b=c 
FIG. 6: The same as in Fig. but for m N = 1 TeV. 



larger than the Z-boson mediated graphs in /i — > e conversion in Gold, for heavy neutrino masses tojv ~ 1 TeV. We 
have checked that for mjv ~ 1 TeV, box graphs due to heavy neutrinos also dominate the process of [i — > e conversion 
in Titanium (not explicitly shown in Fig. Ill|) . Finally, the two lower pannels show the individual contributions 
due to heavy neutrinos A^x,2,3 [blue (solid) line], sneutrinos Nip,... ,12 [red (dashed) line] and soft SUSY-breaking LFV 
terms [violet (dotted) line]. From these two lower pannels, it is obvious that for heavy neutrino masses tojv ;J 1 TeV, the 
soft SUSY-breaking effects dominate the CLFV formfactors, which are tagged with the superscripts SB in Appendix [Cl 
Instead, for m^r < 1 TeV, heavy- neutrino effects start becoming the leading contribution to the CLFV observables 
associated with the muon. Note that the green (horizontal) lines in the middle and lower pannels serve as reference 
values obtained by assuming that only the G^' R formfactors contribute to the amplitudes. 

An important consistency check for our numerical analysis has been to analytically show that all soft SUSY-breaking 
effects on the formfactors (|C.4j) . (|C.8|) and (|C.14j) vanish in the limit of degenerate charged slepton masses. On the 
other hand, RG effects from Mgut to Mz induce sizeable deviations to the charged slepton mass matrix from the 
unit matrix. As a consequence, unitarity cancellations due to the so-called GIM mechanism become less effective in 
this case and so render the SB part of the formfactors, such as F^ B and F^ B , rather large. 

Another essential check was to show that under the assumptions adopted in [ZiJ , our formfactor F^(^ given in (|C.7j) 
reduces to ^2W-P^ } where F£ is one of the formfactors defined in (6) of 44], which in turn can be shown to vanish. 
The assumptions in [44[ are: (i) the standard seesaw mechanism with ultra-heavy right neutrinos, (ii) no charged wino 
or higgsino mixing, and (iii) the dominance of the wino contribution. Under these three assumptions, the interaction 

vertices occurring in the formfactor F^^ simplify as follows: 

AmA> Bhi?A ~^ ~ u ik, C\ B , C\ B , C\ B , C\ B --5 kk > , V* k R -> c 2 w , (4.6) 

where A, B assume now the restricted range of values k,k' = 1,2,3 and U is a 3 x 3 unitary matrix. Given the 
simplifications in (|4.6[) , we recover the expression of [4J], resulting in the replacement: F^{z ~* ^ip-^L- The above 
non-trivial checks provide firm support for the correctness of our analytical and numerical results that we have 
presented in this section. 
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FIG. 7: Numerical estimates of B(fi —¥ ey) [blue (solid)], B(fj> —¥ eee) [red (dashed)], R^l [violet (dotted)] and R^ [green 
(dash-dotted)], as functions of B(/j, —t &y) (left pannels) and the heavy neutrino mass scale mjv (right pannels). In all pannels, 
the Yukawa parameter a was kept fixed by the condition B(ji — > ej) = 1CP 12 for tun = 400 GeV, and tan j3 = 10 was used. 



The upper pannels display numerical values for the Yukawa texture (|2.4|) . with a — b and c — 0, and the lower pannels for the 
Yukawa texture (12.51). with a = b = c. 




m N /GeW m N /GeV 

FIG. 8: Contours of the Yukawa parameters (a,b,c) versus mjv, for B(fi — > 67) [blue (solid)], B(/x — >• eee) [red (dashed)], 
[violet (dotted)] and R^ [green (dash-dotted)], where a and mjv are determined by the condition B(fj, — > ej) — 10 _1 . 
All contours are evaluated with tan ft — 10 and for different Yukawa textures, as indicated by the vertical axes labels. 
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FIG. 9: Contours of the Yukawa parameters (a, b, c) versus tun, for B(t — > cy) [blue (solid)], where tan /3 = 10 and a and tun 
are determined by the condition B(r — s- 67) = 1CP 9 . No solutions have been found for B(t — > eee) and B(r — > e/i/i). 
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FIG. 10: Numerical estimates of _B(p — > e-y) [blue (solid)], B(fj, — > eee) [red (dashed)], i?^, [violet (dotted)] and ii^" [green 
(dash-dotted)], as functions of tan j3. The upper pannels are obtained for mjv = 400 GeV and the lower pannels for mjv = 1 TeV. 
The left pannels use the Yukawa texture (|2.4|) . with a — b and c = 0, and the right pannels the Yukawa texture (12. 5|) . with 
a = b — c. In all pannels, the Yukawa parameter a is determined by the condition B(fi — > e7) = 10 -12 . 



V. CONCLUSIONS 



We have analyzed Charged Lepton Flavour Violation in the MSSM extended by low-scale singlet heavy neutrinos, 
paying special attention to the individual loop contributions due to the heavy neutrinos ^1,2,3, sneutrinos -/Vi ) 2,...,i2 
and soft SUSY-breaking terms. In our analysis, we have included for the first time the complete set of box diagrams, 
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FIG. 11: Numerical estimates of the ratios J?2(m — > eee) 



R£ l and Rs u , 



as functions of mjv- The Yukawa parameter a is 



fixed by the condition B(fi 
complete evaluation of i?2(^ 

lines are evaluated keeping only the magnetic dipole formfactors G 7 and G". The two middle pannels provide a formfactor 



erf) = 10* , for ttin = 400 GeV and tan/3 = 10. In the upper pannel, thick lines give the 
-5> eee) [blue (solid)], Rj 1 [red (dashed)] and R$ u [violet (dotted)], while the respective thin 



analysis of it^A* 



and -R^ u , in terms of contributions due to G 7 and F-y [blue (solid)], Fz [red (dashed)] and box 



formfactors [violet (dotted)]. The lower two pannels show the separate contributions due to heavy neutrinos N [blue (solid)], 
sneutrinos N [red (dashed)] and soft SUSY-breaking LFV terms [violet (dotted)]. The green (horizontal) lines in the middle 
and lower pannels give the predicted values obtained by assuming that only the G^' R formfactors contribute to the amplitudes. 
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in addition to the photon and the Z-boson mediated interactions. We have also derived the complete set of chiral 
amplitudes and their associate formfactors related to the neutrinoless three-body CLFV decays of the muon and tau, 
such as /j, — ¥ eee, r — > fififi, r — > e/i/i and r —¥ eeji, and to the coherent [i — ¥ e conversion in nuclei. Our analytical 
results are general and can be applied to most of the New Physics models with CLFV. In this context, we emphasize 
that our systematic analysis has revealed the existence of two new box formfactors, which have not been considered 
before in the existing literature of New Physics theories with CLFV. 

Our detailed study has shown that the soft SUSY-breaking effects in the Z-boson-mediated graphs dominate the 
CLFV observables, for appreciable regions of the i/rMSSM parameter space in mSUGRA. Nevertheless, there is a 
significant portion of parameter space for heavy neutrino masses win ~ 1 TeV, where box diagrams involving heavy 
neutrinos in the loop can be comparable to, or even larger than the corresponding Z-boson-exchange diagrams in 
fi — > eee and in \x — > e conversion in nuclei (cf. Fig. 1111) . In the same kinematic regime, due to accidental cancellations, 
we have also observed a suppression of the branching ratios for the photonic CLFV decays /i — > e*y, as well as for 
the decays r — ¥ e^j and r — > /17. As mentioned in the Introduction, such a suppression in low-scale seesaw models 
is a consequence of a cancellation between particle and sparticle contributions due to the approximate realization of 
the SUSY no-go theorem due to Ferrara and Remiddi [43|. Instead, in high-scale seesaw models such cancellations 
can only occur for particular choices of the neutrino- Yukawa and Majorana-mass textures [39l l52j. Hence, the results 
obtained within supersymmetric low-scale seesaw typc-I models, with mjv .~ 10 TeV, corroborate the original findings 
in |34[, where the usual par adig m with the photon dipole-moment operators dominating the CLFV observables in 
high-scale seesaw models [37l |38| gets radically modified, such that — > eee and /1 — > e conversion may also represent 
sensitive probes of CLFV. 

We have found that unlike heavy neutrinos, CLFV effects induced by sneutrinos remain subdominant for the 
entire region of the mSUGRA parameter space. In addition, the perturbativity constraint on the neutrino Yukawa 
couplings h. v up to the GUT scale renders the quartic coupling contributions of order (h y ) 4 small. The present study 
has focused on providing numerical predictions for relatively small and intermediate values of tan/3, i.e. tan/3 < 20, 
where neutral Higgs-mediated interactions constrained by the recent LHCb observation of the decay B s — > n[i are not 
expected to give sizeable contributions. A global analysis that includes large tan /3 effects on CLFV observables and 
LHC constraints will be given in a forthcoming communication. 
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Appendix A: INTERACTION VERTICES 

In this appendix, we list the Lagrangians describing the interaction vertices required to calculate the transition 
amplitudes for the CLFV processes under study. The corresponding interaction vertices for the SM and the MSSM 
are obtained by adopting the conventions of the public code FeynArts-3 . 3 . FVMSSM . mod. The Lagrangians of interest 
to us include: 

1. Vertices from 2HDM sector of the MSSM involving SM particles only, 

^ U H-+ H - C - = ^^V* j d j (t m d ,P L + t0 1 m Ui P R )u i H-+ll. C . (A.l) 

Here H~ is the negatively charged Higgs scalar, V is the Cabibbo-Kobayashi-Maskawa matrix, and m Ui are the 
quark masses and c w — cos 8 W . 

2. Vertices of singlet neutrinos in the v R SM sector of the MSSM, 

£enG-+H.c. = Jt^ B ia ei( -m e% P L +m na P R jn a G~ + H.c. , (A. 2) 

£enw- + H.c. = ~^B la ea^PLUaW- + H.c. , (A.3) 

Cnnz = -^C ab n a ^P L n b Z^ . (A.4) 

Here n a and m na denote the neutrino mass-eigenstates and their respective masses and B and C are lepton flavour 
mixing matrices defined in [2^| and [Hj]. The matrices B and C satisfy the following set of relations: 

Bl a B\, a = Sw, C ac Cb c — C a b, BibCba — Bi a , B\ a Bib = C a b, 

m a CacC b c = 0, m a BibCba = °> m a Bi a Bi> a = . (A. 5) 



3. Vertices from the 2HDM sector of the MSSM involving Majorana neutrinos, 

C-enH- +H.c. = -^^B ia ei(tpm ei PL +tp 1 m na Pii)n a H~ +H.c. (A. 6) 

4. MSSM vertices with sparticles, 

« + H.c. = 9 W di (Vr£P L + V^Pr) x-u a + H.c. , (A.7) 

J + H.c = g w u 3 (v+ m u a L P L + V +Z L Pr) xLda + H.c. , (A.8) 

A + C — k+a = e Xm7 ,I X m ^-eXm7''xi^> ( A -9) 

£ -x-z + c jfx+z = ^^r(K k L PL + v*; k R Pu)r k z, 

- ^l^V^Pn + V^pAxtZ,, (A.10) 
Cw ^ ' 

c *°z>z = ^xLirfPivS? +rp R v£ R ) x ° k z, (A.ii) 

Ce*ez = g w V^ b e* a id"e b Z„ , (A.12) 

C^o- e + C ¥eS , = g w e 3 (P L V^ L a + P R V^)x°e a + H.c. , (A.13) 

Cu X °u + %, a . = g w Uj (PlV^ + P R V^)x°Ua + H.c. , (A.14) 

+ = 9»di{PLV$% + P R V™ R )x°d a + H.c. , (A.15) 



where 



v jma — fK 77 lA m2 v ij\ IX L)ai J 
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m d . 
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U m 2Vji(R R )* ai 
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w ' 



V** = V* x V kl + \v* m2 V k2 - 6 mk s 
— 7(^3^3 — z m4 zi 4 ) , 
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w ' 



v ab 



K H L)ai{n L ) bi (ttR)ai(n,R)b 

Cyj C-w 



-V2t w Z* ml {R 
1 



e \* 
RJaj 



(m e )j 



\/2cw 

2y/2 



y/2cpM w 
(cwZ m 2 + Sw Z m i){R e L ) 



e \* 
aj 



(™-e)j 

V2cf 3 M w 



-Z m z(R R )* 



"j 



{ m u)j y* (r/u\* 

' /j mA\ sx Llaj 1 



2cw 
V2 



"' y/2s p M w 



(c\vZ m 2 + 7 Sw Z m i)(R l l)*j j- — Z mi (R R 



twZ ml (R R/ 



(m d )j 



"j 



-Zm?,{R\ 



1 W 



LI aj 5 



— -{c w Z m2 - -^s w Z ml )(Rf)* aj - J^^m Z mZ (R^^ 



Rlaj ) 



(A.16) 

(A.17) 

(A.18) 
(A.19) 
(A.20) 
(A.21) 
(A.22) 

(A.23) 

(A.24) 

(A.25) 

(A.26) 

(A.27) 



and c 2w — cos 26 w . The unitary matrices diagonalizing the chargino mass matrix 11 and V and the unitary matrix 
diagonalizing the neutralino mass matrix Z are taken from [69| . The matrices 



R 



J'L 



u ia u ik ' 



R 



fR 



Tjf TjIr* 
u i+3a u ik ' 



(A.28) 



with / = d, u, e, a — 1, 2, . . . , 6 and i, k = 1, 2, 3, quantify the disalignment between fermions and sfermions. Here 
jjih ^ jj/r an( j jjf are un itary matrices that diagonalize the fermion and sfermion mass matrices, respectively. 



5. Sneutrino vertices in the i^rMSSM, 

^eic-N + H.c. = g h 



where 



^NNZ 



9w N*Jl { P ^™ l Mw BlmA + PRBt\)xm + H.C. 

cw 

^C 2 AB N^N B Z, 
cw 

^LC\ B N A ^N B Z^ 
c w 

^C\ B N A id^N* B Z^, 

Cw 



(A.29) 



(A.30) 
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C = -UiA U u L *V m l + -J=-^f—Vrn2Ui +9A Ui^ Za Bla, 

AmA = ^i+eA^^ml + ^2 S pM w Vm2 ^+ 3A U i+3a B la , 

^7 7^*7 7^ 

^7yp*7yl>* 

b AB — 2 iA i+6B> 
~3 1 i> £> 

Cab = ~-jUi+sAUiBi 

Cab = -\uf +6A UtUB • (A.31) 

In the above, U v is the unitary matrix diagonalizing the sneutrino mass matrix. 

Notice that we have factored out the weak coupling constant g w from all interaction vertices defined above. To better 
identify chirality-flip mass effects in the CLFV amplitudes, we have also pulled out a factor mi/ 1 (\/2cpMw) from the 
interaction vertex Bf mA . 



Appendix B: LOOP FUNCTIONS 



The CLFV amplitudes are expressed in terms of leading-order one-loop functions. We expand the loop functions 
with respect to the momenta and masses of the external charged leptons, while keeping only the leading non-zero 
terms. The leading terms may then be expressed, in terms of the dimensionless loop integrals 



J ni n 2 . .. n k (^l>^2, ■■■ Ak) — 



^f-D/2 , d D £ (flyn 



{M^)- D ' 2 - m+ T.^ J {^) D n- =1 (^ 2 - m?)«* 

i(-l) m -Ei n ' / n 2 X2-D/2 foo dxx D/2-l+m 



^m- w ) J n fe -i (x + \i) ni ' { ■ } 



where are loop particle masses, i%i are the exponents of the propagator denominators, A^ = mf/M^ are dimension- 
less mass parameters and n is 't Hooft's renormalization mass scale. We choose u to be M w , even though any other 
scale can be chosen equally well as a reference scale for any of the integrals. For the amplitudes we have been calcu- 
lating, the integrals are either divergent and satisfy m + 2 — JV rii = 0, or they are convergent with in + 2 — £V rii < 0. 
For convergent integrals, we may set D — 4, whilst for divergent integrals we take D = 4 — 2e. We pull out a factor 
i/(47r) 2 from all integrals. Thus, we have for finite integrals, 

^n 1 n 2 ...n Ai (Al,A2,...,A/ s ) EE 2 J™ lTl2 ... nk (Ai , A2, . . . , Afe) . (B.2) 

Instead, the divergent integrals are written down as a sum of a divergent+constant term and a finite mass-dependent 
term: 

7T 1 „ 2 ..., lfc (Ai,A2,...Afe) ee (i+const + J™„ 2 „ fc (A 1; A 2 , . . . , A fe )) . (B.3) 

In the CLFV amplitudes, the divergent+constant terms vanish in the total sum, or as a result of a GIM-like mechanism. 
Therefore, all CLFV amplitudes can be expressed in terms of finite mass dependent functions J™ lTl2 (Ai,A2,...), 
which we call them the basic integrals. The CLFV amplitudes and the corresponding formfactors considered here are 
described by 9 basic integrals, four for the photonic amplitude: J^, Jg 1; J\ x and J| 1; three for the Z-boson amplitude: 
J n , J nl and Jin, and two for box amplitudes: J nll and J nil . 
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Appendix C: ONE-LOOP FORMFACTORS 



Here we present the complete analytical form of the CLFV formfactors F 1: Fz and -Fbox defined in Section HlT| in 
the Feynman-'t Hooft gauge. In the following, the usual summation convention over repeated indices is implied. The 
interaction vertices and loop functions used here are given in Appendices [X] and |B] respectively. 



1. Photon Formfactors 



The formfactors i*]r, F^, G 7 and G 7 may be explicitly written as follows: 



where 
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yO£L yOiL* 
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.Z-Boson Formfactors 



The formfactors F z and may be decomposed as follows: 
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3. Leptonic Box Formfactors 



The leptonic box amplitudes are expressed in terms of the chiral structures: W^l hT^l^ [cf. (|3.8j . There are 
two distinct contributions to the chiral amplitudes. The first one has direct relevance to the original structure given 
above and we denote it with a subscript D. The second contribution comes from a chiral amplitude of the form 
hT^l ?r^2 , which contributes to the original amplitude I'T^l hT^l^ , after performing a Fierz transformation. 
This Fierz-transformed contribution is indicated with a subscript F. More explicitly, the leptonic box formfactors are 
given by 



R LL 


nLL 

— n ev,D 


+ Bey^p, 




r>RR 
B £V 


nRR 

D ev,D 


i nRR 

' n ev,Fi 




rjLR 
B IV 


r>LR 

— D iV,D 


1 nLR 

~ 2 U tS,F> 




nRL 


nRL 

— n ev,D 


1 Rfl-L 
~ 2 n lS,Fi 




nLL 
a tS 


nLL 

n es,D 


+ 2~ B lkF "+ 


rLL 


r>RR 
U £S 


nRR 

— a £S,D 


+ 2 is,f "+ 


" 2~ n eT,Fi 


tdLR 
B £S 


nLR 

n es,D 


La £V,Fi 




nRL 
B £S 


nRL 

n es,D 


o nRL 
La £V,Fi 




f>LL 


nLL 

— D IT,D 


1 nLL , 
~~ 2 °IT,F i 




tdRR 


nRR 

£T,D 


1 R il.R 
- 2 -°€T,_F " 





(C.9) 



The direct and Fierz-transformed contributions to the formfactors are related by the exchange of outgoing leptons 
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/dL,1 rA,1* i r£. 2 R fi,2* unSJ R 
\ D l 1 kB D l' 2 mB J - > l 1 kB JD l 2 mB)\ JD l'mA r> i 



B V rr, A^lhA )J\\l\^Xk ! Ax m 1 A^ , Ajy ), 



£,L,1* 



/'mA Ik A 
jR,2 R L,2* 



mi' mi 

9 r 2 If 2 



•^fiii (Axfc j A^ m , A^y , A^ ), 



/ R iJ,l f>L,l* f,R,2 R L,2* w„fi,l R 
y D l 1 kB D HmB " r D l\kB D limBly V mA 13 . 



Ik A 



nR,2 R L,2*^ 



vA^A^m "^liii (A** ! Ax m , A^ A , A^ b ), 



/ R H,l R -R,l* , R R,2 R il,2* wni,l R L,1* R L,2 R L,2*s, 
\ D l 1 kB D l. 2 rnB^ D l 1 kB D l. 2 m,BI\ c 'Vm,A D lkA ^ D VmA D lkA > 
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2 ^2 \/ A Xk A X m Jllll(. A Xki A X m i A N A i A N B )i (CAS) 
C /3 1V1 W 



r>LL,SB _ frOmfrOtR*frOtRfrOtR* rl I \ \ \ \ \ 
n tV,D ~ ^ v l 1 mb v lma v l'na v l 2 nb ^llll x A X°,. ' A x£. ' A e» > *e b ) 



c,fMRfMR*fMRfMR* R x tO (\ x x x \ 
_ zv l 2 mb V lma v l'na V hnb y A X° m A X° k \ A X° m ' A X° n ' A e« ) A e 6 ), 

U RR,SB _ fML -rML*-rML-?ML* jl / x x x x x 
IV, D — - V hmb V lma V l'na V l 2 nb J llll \ A X ( L ' A x£ ' A e. ' A ^b J 

c,irOlL frOlL*frOLL-rML* R x 7 /x x x \ x 
ZV l 2 mb V lma v l'na V hnb \J A X° m A X° k J llll \ A X° m ' A X° n ' A e« ' A e (J> 

D Lfl,SB _ o f>0£L f>0£il*f>0£ilf>0£L* /\ x 7 / x x x x x 
^V.D — ^lim^lroa H'no H 2 n6 y *m A X? J im lx ' A x£ ' Ag a ) A e J 

, f/Ott f/0«iJ*i>0Mf/0£L* il fx x x x \ 
"•" v l 2 mb V lma "i'no v hnb J llll v A X^ ' A X° ' Ag a i A et J: 



^V.D - 21/, 



J IS,D 



roeR-rML*-r r oiL-r r MR* R x T o /x x x x \ 

hmb V lma V l'na V l 2 nb V/ A xSk A x£ J llll V A xSk ' A x£ ' A e. ' A e J 

, f>0£iJf>0£L*f>0£Lf>0«fl* rl /x x x x \ 
' v l 2 mb v lma v l'na V hnb 1111 \ Xm ' Xn 1 ^ a ' ^6 / ' 

rOtL ■Cr0tR*fr0tL-?r0lR* R x tO fx x x x X 

hmb v lma v l'na v l 2 nb U A x° m A X k 1111 v A X™ ' A X° ' Ag a ) A e b ) 

ofrOiRirOlRr-TrOiL-rrOeL* R x tO /x x x xx 

_ ZV l 2 mb V lma V l'na V hnb \J A X° m A X k HH ^ A X° m ' A x1 ' Ag a ) A £ J ) 

R -RR,SB _ 0T /(«-R fML*fMRirOtL* R x tO i x x x x x 
^5,13 _ - ZV hmb V hna V l'na V l 2 nb y A Xm xj lllll ft ' X° ' ^ a ' ^ 6 ' 

O fr0£L frOiLtfrMRfrOeR* R T tO x x x XX 

~ ^^mfc^ma V l'na V hnb \J A X° m A X° k J HH v A xS« ' A X^ ' Ag a ) A e J: 

R ifl,SB _ O T>0£ijT/0£il*f>0£Lf>0tt* rl /x x x x x 
£S,D — lV hmb V lma V l'na V l 2 nb J llll l A X^ ' A X° n ' A e. ' A ^b ) 

j_ of>0££ fMR*iML-TrOlR* jl /x x x x x 
+ Z|/ l 2 mFIma "i'na ^iint J llll l A X^ ' A X^ ' A e. ' A e J J 

pfli,SB _ csfML fML*-TMRfMR* T l /x x x x x 
n iS,D " zv hmb v lma v l'na v l 2 nb J llll x A X„ ' A x£ ' A e. ' A e b ) 



l OT/M" T/WLtT/WKi/U«L* rl / x x x x x 

"+" Z H 2 m/i r ima M'na ^int J llll l A xSi ' A X^ ' A e. ' A Sb J: 

n LL, SB _ O t>0(!L f/0^fi*f>0^Lf/0^fi* /x x tO ^x x x x x 
n iT,D — / nim6 )/ ima v l'na V l 2 nb \l A X° m A X k 1111 v A X^, ' A X° ' Ag a J A e b J 

I OT>0^T>0£il*f>07?if>0tt* /T ^ tO /X X X_ XX 

' zv l 2 mb v lma v l'na v hnb V/ A x° m A X k HH V A xSk ' A x" ' Ag a ; A e J j 

n itit,SB _ O T>0(!fl T>0£L*T>0£flf>0(!L* /T T 7 /x x x x x 
^T.n ~ - ZV l im b V lma V l'na V l 2 nb \J A ft A X° k J 11U \ A X° m ' A X^ ' Ag a ) A e J 

+ 2Vl° 2 mbVlma *S ^ A x£> A X? J llll ( A X^ ' A X^ ' A ^a ) A e 6 ) • (C14) 

4. Semileptonic Box Formfactors 

Semileptonic form factors have only direct contributions, with the following N, N and SB content: 

r>LL r/ LL > N i r>LL,N . U LL, SB 

B dV ~ B dV + B dV + B dV 

r>LL tjLL,N U LL,N R LL, SB (n -, rx 



and 



R JCY _ R XY : SB pXY _ R XF,SB /p 1 A x 



for (X, y, A) 7^ (L, L, V). The N and iV contributions are given by 



1 H — a — - ) J2n(l, A„ a , A U!) ) + 2A„ a A U(J J2 11 (l, A„ a , A U() ) 



+ 2r2" A " a ^"^ini^' A ^+' A ™a> A u b ) — 2£2 A n a ^u b Jiin(^, A ff +, A„ a , A„ 
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B 



LL.N 
uV 



LL,N 
dV 

LL,N 



uV 



4t* 



A« a A Ui) J 211 (A^j+ , A„ a , A Ub ) 



S /'a- B ra(^ / *)d2b(^ / )rflfc 



^4-1 — ""^ — ^J2ii(l! A Wo , Ad,,) - 2A na Ad b J2ii(l; A„ a , AdJ 



+ 2^™ a ^ <ii '^i 111 ^' AH+ ' ^«a7 Ad b ) — 2^™ ti ^ di ''^i 111 ^' AH+ ' Adf >) 

+ ^A na Ad b J2h(Ah+) A„ a! Ad 6 ) 

— _ ^T/-dflT>-dR*/T7-^,l*T/-^,l I T/-^,2* T >-^i?,2x t1 fx \ X \ \ 

~~ 9 r tiifca V d 2 ma \ V IkA V l'mA ~ tV lkA V l'mA I J llll \ A *x h ' A X™ I A N A I A »„ 7 



(C.17) 



- -KifeaKam^C^feA ' t'mi' + ^ZfcA ' ^'mA ) V 7 A Xfc A X™ J llll ( A Aj fc , A Sm , A^ A , Aj^ ) . (C.18) 



XY SB 

leptonic formfactors ' , by making the replacements: £ — > d, l\ — > d, I2 — > d, e — > d and £ — > u, l\ — > u, I2 — > it, 



The SB form factors and -B,fJ' SB , with X = L, R, Y = L, R and A = V, S, T, are obtained from the direct 



e — y Uj in both the interaction vertices and the arguments of the J- loop functions that carry the index b in (|C.14[) . 
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